Momentum— space engineering of gaseous Bose-Einstein condensates 
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We show how the momentum distribution of gaseous Bose-Einstein condensates can be shaped 
by applying a sequence of standing-wave laser pulses. We present a theory, whose validity for was 
demonstrated in an earlier experiment [L. Deng, et al., Phys. Rev. Lett. 83, 5407 (1999)], of the 
effect of a two-pulse sequence on the condensate wavefunction in momentum space. We generalize 
the previous result to the case of N pulses of arbitrary intensity separated by arbitrary intervals and 
show how these parameters can be engineered to produce a desired final momentum distribution. 
We find that several momentum distributions, important in atom-interferometry applications, can 
be engineered with high fidelity with two or three pulses. 

PACS numbers: 03.75. Gg,67.85.Hj,03.67.Dg 



I. INTRODUCTION 

The ability to create gaseous Bose-Einstein conden- 
sates in the laboratory [l|-[5[ has led to new vistas in the 
field of atom interferomctry. This is particularly true 
when laser light is used to manipulate atoms to produce 
matter-wave interference patterns. Atom interferome- 
ters using light gratings acting on matter waves have 
been used in a variety of fundamental studies such as 
how large a composite object can display interference ef- 
fects [H , decoherence studies @, HI , origins of phase shifts 
under various circumstances, properties of Bose-Einstein 
condensates [ill-fill], and testing the charge neutrality of 
atoms [l2j . Atom interferometers are also at the heart of 
a host of practical devices used for making precision mea- 
surements. These include gravimeters, gyroscopes, and 
gradiometers all of which have important applications in 
precision navigation [T§ - [l5j ]. Such interferometers also 
have applications in atomic physics such as atomic po- 
larizability measurements and Casimir-Poldcr potentials 
for atoms near surfaces [p6| . More in-depth information 
about the uses of atom interferomctry can be found in 
Ref. [13 ■ 

One of the crucial elements of an atom interferom- 
eter is initial state selection of the atoms and these 
states are generally states of localized momentum fl7j . 
Momentum-state selection techniques are quite varied 
and can range from using a pair of collimating slits that 
select thermal atoms with limited transverse momentum 
to Bose-Einstein condensation of atom clouds via tech- 
niques that are a combination of laser cooling and trap- 
ping and evaporative cooling. In this paper we show 
how applying a sequence of short-time, standing-wave 
laser pulses to a Bose-Einstein condensate (BEC) can 
be used as a tool for the state selection step of a atom 
interferometric experiment. We will also see that prepa- 
ration of several important classes of momentum states 
can be achieved through the application of just two or 
three pulses. 



Previously a sequence of two short-time, standing- 
wave pulses (sometimes called Kapitza-Dirac pulses) has 
been used as the beam splitter in the experimental re- 
alization of a Michelson atom interferometer for a BEC 
formed on an atom chip [Lsl ]. Optimization of the two- 
pulse sequence was determined by studying a two-state 
truncation of the Raman-Nath equations [lj|. Other 
related beam-splitter techniques which produce specific 
momentum orders [20U22I I have been studied. However, 
these techniques represent a different strateg y [l| in the 
standing-wave control of atomic motion. 



This paper is organized as follows. In Section [TU we 
derive a theory for how the condensate wave function is 
changed following the application of two pulses of differ- 
ent intensity and separated by a time interval on the or- 
der of the Talbot time. The Talbot time is Tt = h/E reco u 
where E reco u is the recoil energy of the atom for the laser 
light used to make the standing-wave pulses. This theory 
is composed of two parts, (1) the effect of a pulse on the 
condensate wave function and (2) evolution of the wave 
function between pulses. These two elements are used 
to follow the steps of a two-pulse sequence to derive an 
expression for the amplitude for atoms to jump from the 
zero-momentum initial state to an arbitrary momentum 
order. In Section IIII1 we derive two general symmetries 
of the two-pulse amplitude and some values for special 
time interval values. This section also presents a physi- 
cal interpretation of the amplitude. This interpretation 
is especially useful for generalizing the two-pulse result. 
Section IIVI contains a derivation of the general TV-pulse 
amplitude. Section [V] presents the least-squares method 
for designing a general sequence of N pulses which gives 
a specified momentum probability distribution. Exam- 
ples of two-pulse and three-pulse sequences are given for 
some important momentum distributions. Conclusions 
are presented in Section IVI1 
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II. TWO-PULSE THEORY 

In this section we derive how the condensate wave 
function is changed when two short-time, standing-wave 
laser pulses arc applied. We assume that each pulse is 
square-shaped in time and that its duration, St, is short 
enough such that there is no appreciable spontaneous 
emission during the pulse and short enough so that are 
no effects of atom-atom interactions during the pulse. 
This will happen if 5t < K/T where T is the natural line 
width of the excited state. We also assume that the entire 
pulse sequence is short compared to Ti/ '/j,, where \i is the 
chemical potential, i.e., short enough so that there are no 
effects due to the interaction during the pulse sequence. 
It is important to note that, under these assumptions, 
many-body effects can be neglected. 

We assume that the first pulse is applied at time t\ 
and has duration 8t\ and that the second pulse is ap- 
plied at time ti and has duration St2- We also assume 
that the time interval, £2 — ii, is small enough so that 
the maximum distance that condensate atoms in non- 
zero momentum states move after the first pulse is small 
compared to their de Broglie wavelengths. This is the 
Raman-Nath regime. 

Our approach will be to use these conditions to approx- 
imate the effect that a single pulse has on the conden- 
sate wave function and separately to approximate how 
the wave function evolves between pulses. With these 
effects in hand, we can then follow the steps of the pulse 
sequence to determine the overall effect of the full pulse 
sequence on the wave function. 

A. Effect of the first pulse 

Consider atoms in a Bose-Einstein condensate inter- 
acting with a pair of linearly polarized, counterpropagat- 
ing laser pulses. Each atom is modeled as a two-level 
system having a ground state \g) and an excited state 
\e). We denote the position of the atom's center-of-mass 
(CM) relative to an arbitrary coordinate system by r a 
and the position of the atomic electron relative to the 
CM by r e so that the position vector of the electron in 
the arbitrary system is r Q + r e . 

The Hamiltonian for a single atom interacting with the 
light and including the CM motion is given by 

H = H^(r a )\g)(g\+H^(r a )\e)(e\ 

+ E g \g)(g\ + E e \e)(e\ + V lascr (r OJ i) . (I ) 

Where 

H^(r a ) = ^ + V^ p (r a ), k = g,e (2) 

are the energies associated with the CM motion of the 
atom in the ground and excited states. The difference in 
trap potentials derives from the different magnetic mo- 
ments for the two internal states. Both potentials are 
assumed to be harmonic here. 



The term Vi aser is the usual dipole interaction and is 
written as 

Mascr (l"a, 

t) = 2hfl cos (k L ■ r a ) f (*) cos (wii) 

x(|e>fo| + |ff><e|), (3) 

where the standing-wave laser field is assumed to have 
the form 

E (r a ,t) = 2£ Q f (t) e L cos (k L • r tt ) cos (ui L t) . (4) 

The laser-field amplitude, frequency, and wavevector are 
denoted by £q, lol, and k^, respectively, while is the 
laser polarization vector. The laser-field amplitude enve- 
lope, / (t), is assumed to be a rectangular pulse centered 
at an arbitrary time t% with width 8t\. The factor f2o in 
Eq. ([3]) is the single-photon Rabi frequency given by 

hn = e£ (e\r e ■ e L \g), (5) 

where e is the electron charge. 

The condensate orbital is represented by a spinor wave 
function of the form 

= ip g {r a ,i)\g) +ip e (r a ,t)\e) . (6) 

These components satisfy the multicomponent Gross- 
Pitaevskii (GP) equations 

ih^ = H^ip g (r a ,t)+E g ij g (v a ,t) 

+ 2hil cos (k L • r ) / (t) cos (ui L t) tp e (r , t) 
+ gN (|^ g (r Q ,t)| 2 + \i/> e (r a ,t)| 2 ) ^ g (r OJ t) , 

(7) 

and 

ih^- = H^ e) ^ e {r a ,t) + E e ip e (r a ,t) 

+ 2hn cos (k L ■ r tt ) / (t) cos (uj L t) ip g (r a ,t) 
+ ,giv(|V g (r a ,i)| 2 + |^e(r a ,i)| 2 ) 1>e(r a ,t). 

(8) 

The solution of these equations, over the time of the 
laser pulse can be approximated if (I) the pulse time 
is short enough so that 6ti <C h/fi where fi is the con- 
densate chemical potential, and (2) if the single-photon 
Rabi frequency is small compared to the detuning from 
resonance, i.e., f2o *C A. The detuning is defined by 
fiA = E e — E g — %uii,. We will see that the solution is 
valid even for a strong pulse. 

By approximately solving the multi-component GP 
equations over the duration of the pulse, the details of 
which arc given in Appcndix[Xl the effect of a short-time, 
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strong-field, standing-wave laser pulse on the condensate 
wave function is given by 



(r ,ii+) « e 



xe 



ifi2<5ti cos(2k^'r n ) 



where 



n - n§ 

" 2 =2A 



•'•^ fl (r a ,fi_). (9) 



(10) 



is the two-photon Rabi frequency. 

It is possible to represent the effect of the pulse in 
momentum space by using the Bessel generating func- 
tion m. 



(11) 



n— — oo 



Letting i = je - ™ CI, ' ro anc l z = f2 2 <5ii gives 



itizSt! c OS (2k L -r a ) = i n J n (n 2 St l )e 2nikL - r '. (12) 



n— — oo 



Hence, the wave function just after the pulse can be writ- 
ten as 



(r a ,ti+) « e 



J2 i n Jn{Sl25h)< 



,2nikr r c 



71— — OC 



The 



„2mk.L -r a 



(r 0) *i-) ■ 
(13) 



(r a ,ti_) factor in the n term in the 



above sum is the wave function of the original zero- 
momentum condensate kicked into a momentum state 
centered at p„ = 2nKk.L and expressed in the position- 
space representation. Thus the amplitude for a conden- 
sate atom starting from zero momentum and kicked into 
momentum state 2nKk.L by the pulse is i n J n (f2 2 (5ii). We 
now turn to the evolution of the condensate wave func- 
tion between pulses. 



B. Evolution between pulses 

Between the first laser pulse at t = t\ and the second 
at t = t 2 , the evolution of an atom in momentum state 
p„ = 2nKk.L can be approximated as a free particle under 
conditions described below. The phase of such an atom 
thus evolves as e~ lEn ^ t ' 1 ~ tl ^ h where E n = p 2 l /2m is its 
kinetic energy and the atom moves with velocity v n = 
2nh\s.L,/m. For non-zero momentum states, as long as 
the number of atoms outcoupled into them is too small to 
be detected experimentally or the time scale over which 
they evolve is small compared to %/ mean-field effects 
can be ignored. Thus we can write the condensate wave 



function during the time t\ < t < t 2 as 



-ti)/n 



71 — — OO 

x e 2nik ^d> g (r a -v n (t-t 1 ),t 1 ) (14) 

It is clear that, although the above equation implies 
that there is a finite probability for atoms to be outcou- 
pled into a momentum state that is any multiple of 2Kk.L, 
above some maximum order, n max , there will be too few 
atoms present to be detected experimentally. Here we 
will assume that the Raman-Nath approximation holds, 
that is, the only momentum orders appreciably popu- 
lated are ones in which the atoms moved only a small 
fraction of their de Broglie wavelengths during the time 
between pulses. In this case, 



(r 0> t) « e iQ ^ J2 i n Jn(n25h)e- iE - 



(t-ti)/fi 



x e anikl - r -0 i ,(r o> t 1 ). (15) 

Note that we have neglected the motion of the non-zero 
momentum orders over the time interval t — t\. We shall 
assume that the Raman-Nath approximation holds here- 
after. Later, when we consider the iV-pulse case, we will 
assume that this holds for the entire pulse sequence. Next 
we analyze the effect of the second pulse. 



C. Effect of the second pulse 

The effect of the second pulse can be described by ap- 
plying the exponential in Eq. ([§]) to the wave function in 
the Eq. (|15[) . The wave function just after the application 
of the second pulse is 



(r ,t 2 +) ~ e 



ifi 2 (<5ti+<5* 2 ) 



X 



-iE„(t 2 -ti)/S 



x e 2nik L r a e *H 2 8t 2 cos(2k I , ; a ) ^ ^ ^ _ 



(16) 

The exponential describing the effect of the second pulse 
can also be expanded in a truncated series of Bessel func- 
tions and we have. 

cP g (r a ,t 2+ ) » e ^(s tl +st 2 ) 

Umax TL max 

x £ £ i n+n 'j n (n 2 Sh) 

n=—ri m , ax n'=—n max 

x J n , (Q 2 5t 2 )e- iE " ( ' t2 - tl ^ h e 2 ( n+n ') ikL - ra 
x 4> B (T a ,ti). (17) 

Changing summation indexes as m = n + n' we can write 

4> g (r a ,t 2+ ) » e ^(5t 1+ st 2 ) 
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m— — 2n 7 - nax 



^•'•^ (r o ,t0 ; 



(18) 



where 



A m = Z m Jm-n(^2Sh)J n (fi 2 Sh)e- iE ^-^/ H . 

7l= — OO 

(19) 

Note that we have extended the limits of the summation 
back to infinity. This is possible because, in the Raman- 
Nath regime, the values of J n (£l 2 5t) for n > n max are 
sufficiently small that we incur little error in including 
these extra terms. So, finally, we have 



(r a ,i) ~ e 



ifl 2 (St 2 +St 1 ) 



A m e 2mi ^^ g (v a M)- 



m— — oo 



(20) 



The quantity A m is the probability amplitude for an 
atom to be in momentum state p m = 2mKkr 1 . 



III. GENERAL FEATURES OF THE SOLUTION 
FOR TWO PULSES 

Before deriving some general features of the two-pulse 
amplitude, it will be convenient to rewrite Eq. (|19p in a 
form where the time interval is measured in appropriate 
units. To this end we define the period, Tt, of free oscil- 
lation of a free particle whose momentum is pi = 2Kk.L 
such that: 



Pi 



^- T T /h = 2tt. 
2m' 1 



This is the Talbot time defined above. In this case 
E n (t 2 -h) 2 (pi \ (t 2 -h) 



(21) 



2m / h 



2nn 2 (t 2 -t 1 ) _ o „ 2 



2tt/W, (22) 



where we have introduced the time interval parameter 
as Pi = (t 2 — tij/Tx, that is, the time between pulses 
measured in units of the Talbot time. Also defining the 
pulse area parameters as ot\ = ^l 2 8t\ and a 2 = tt 2 5t 2 , we 
can write Am as 



A m (a 1 ,a 2 ,/3 1 ) = i m £ J m _„ (a 2 ) e - 2 -" 2 A J n ( ai ) . 

n — — oo 

(23) 

We will use this form of the two-pulse amplitude to 
demonstrate some its general features. 



A. Symmetry of the momentum distribution 

There are two general symmetries exhibited by the ex- 
pression in Eq. (|2"3"|) . The first symmetry is that the am- 
plitudes for opposite momentum orders are equal for a 
given pulse sequence, that is 



A m (a-L,a 2 ,Pi) = A_ m (ai,a 2 ,/3-i). 



(24) 



To show this symmetry, we write the expression for 
A_ m (ai,a 2 ,/3i) which is 

oo 

A_ m (ai,a 2 ,/3i) = i~ m ^ J-m-n(oi2) 

n— — OO 

x e- Wft J„(ai). (25) 



Using the Bcsscl generating function [23[, it is easy to 
show that J-n(z) = (—) n J n (z) so that 



A_ m (ai,a 2 ,j8i) = i m 



Jm+n (ct2) 



X e 



= I 



-2-irin fii i \n 



(-)" J- n (ai) 
^m+n («2) 



n— — oo 



x e 



-2Trin 2 /3 1 



J-n (ai) • 



Changing the summation index to n' = — n we have 

oo 

A- m (ai,a 2 ,Pi) = i m ^ J m~n' {a 2 ) 

n' — — oo 



x e 



-27ri(-n') 2 /3 



1 J n / («i) , 



and so 



A- m (ai, a2, j8i) = A m (ai, a 2 , . 



(26) 



Thus all distributions are exactly symmetric with respect 
to momentum order under these conditions. Physically 
when atoms make a transition to a non-zero momentum 
state they must absorb photons from one laser beam and 
emit into the other beam. Equation ([26]) holds because, 
for standing-wave laser beams, there is no preference 
for choosing which beam photons are absorbed and into 
which emitted. Absorbing from opposite beams results 
in populating opposite final momentum states. 

The second symmetry is that the amplitude for a given 
pair of pulses separated by an interval f3\ is the complex 
conjugate of the amplitude for an interval of 1 — /?i for 
fixed ax and a 2 . Thus we have 

A m (a 1 ,a 2) l-p 1 )=A* m (a 1 ,a 2 ,p 1 ) ft < 1/2. (27) 

If we only consider interval times /3i < 1/2, it is easy 
to see that this holds in Eq. (H]). Since e -^"^ 2 ^-M = 
e 27 ™" ° l gives the complex conjugate of the exponential 




FIG. 1: (a) (color online) The final distribution of the condensate atoms over the momentum states after the application of two 
standing-wave pulses separated by an interval can be understood in terms of interfering quantum pathways. The amplitude 
for the system to jump from the initial zero-momentum state at the top to the final state at the bottom along the path shown 
is the product of the amplitudes (shown in boxes) for the three legs that compose the path. The total amplitude is coherent 
summation of this composite single-path amplitude over all possible paths, (b) The general case of N pulses is illustrated in 
the right panel. We assume that pulse 1 is applied at ti, pulse 2 at t%, etc. 



in the fi\ amplitude. Given the fact that the Bessel func- 
tions are all real, conjugating the exponential conjugates 
the entire expression. Equation (|27p also implies that 
the probability distribution for two-pulse sequences with 
intervals (3i and 1 — /3i are identical. This result will 
be useful later to constrain the parameter space in the 
momentum space design procedure described below. 



B. Special values of the time interval 

Two important results can be obtained for the quantity 
j4 m (ai, a 2 , Pi) when the interval time has special values. 
The two special cases are (1) (3\ = i, and (2) /3i = 1. We 
consider each of these in turn. 

The first case occurs when the time interval between 
standing-wave pulses equals half of the Talbot time. 
Here, the exponential factor in Eq. (|23|) becomes 

e" 2 "" 2 ' 91 = e-™ 2 = H" 2 = H" . (28) 

The last equality can be seen by noting that the square of 
an even integer is even and the square of an odd integer 
is odd. Thus, the probability amplitude becomes 

A m (a>i,a 2 , - j = i m ^ (-)" J m -n (aa) Jn («i) • 

^ ' 71 — — OO 

(29) 

This expression can be summed exactly as follows. Wc 



write again the Bessel generating function 

OO 

e *«i(*-*) = £ t n J n (ai) 

71— — OO 

OO 

e-WM) = J2 (-f t n ' J* (a 2 ) . (30) 

n' — — oo 

In the second equality, we have let ct\ — > —a 2 and 
used the identity J n (—a) = (— ) J n (a). Equating the 
product of the right -hand-sides with the product of the 
left-hand-side of the above two equations and letting 
m = n + n! yields the following 

e _x (aa _ ai)(M ) = g €»A m L l ,a 2 ^t m 

m— — oo ^ ' 

oo 

= E (-) m Jm(a2-ai)t m , 

m— — oo 

(31) 

where the second equality comes from a direct application 
of the Bessel generating function to the exponential on 
the left-hand-side. Since the equality of the two sums 
must hold for any value of t, the coefficients of t m in the 
sums must be equal. And so, 

A m (a!,a 2 , = i m J m (a 2 - «i) (32) 

One consequence of this is that, for equal-area pulses and 
when the interval between pulses is one-half of the Talbot 
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FIG. 2: (color online) Comparison of the theory presented in the text with the results of the temporal Talbot experiment 
reported in Ref. [llj. (a) The left panel shows the steps of this experiment. Two 100-ns, 1 W/cm 2 , standing-wave laser pulses 
were applied to a 23 Na BEC after allowing the condensate to expand for 1.2 ms. The pulses were separated in time by an 
interval, St, which varied between 0.1 and 11 [is. (b) The middle panel shows the images taken after allowing the condensate 
to expand for a further 6.2 ms. Each column of dots corresponds to a specific final momentum state. The highlighted column 
corresponds to zero momentum while the column just to the right (left) of it corresponds to the 2hk (—2hk) momentum 
state, (c) The right panel shows a comparison of the prediction of Eq. (|19p with the normalized pixel count of the dots of the 
highlighted column. 



time, (that is, when ct\ = a 2 and /3i = 1/2) the ampli- 
tudes for all non-zero momentum states are zero and the 
condensate is unchanged. This effect was verified exper- 
imentally and reported in Ref. . Figure [UJb) shows 
this effect where a double-pulse delay of 5 /j,s results in 
no change to the original condensate. 

In the second case where fix = 1 the exponential in the 
expression for A m (ai, a%, 1) equals unity for all values of 
the summation index n. Thus we can write 

oo 

A m (a 1 ,a 2 ,l) = i m ^ J m -n (cv 2 ) J„ (eti) 

'It — — OO 

= i m J m (a 2 + a x ) (33) 

where the second equality is derived by a method similar 
to that which produced Eq. (|32|) . This result suggests 
that two standing-wave pulses separated in time by one 
Talbot time have the same effect as a single pulse whose 
area is the sum of the areas of the two pulses. 

C. Physical interpretation of the probability 
amplitude 

The general formula for the probability amplitude, Eq. 
([TT)|) . can be understood as the superposition of ampli- 
tudes of multiple pathways from the given initial to the 
given final state. One such pathway is illustrated in Fig. 
[IJa). This figure depicts a particular quantum pathway 
from the fixed initial zero momentum state to a fixed 



final state whose momentum is p m = 2mhk. The ampli- 
tude for an atom to go between these states via the path 
shown is the product of the amplitudes for the three legs 
of the path. The first leg is a momentum jump, caused by 
the first pulse, from the zero-momentum state to the mo- 
mentum state p = p n = 2nhk with amplitude i n J n (a). 
In the second leg of the path, atoms in the momentum 
state p — p n , whose energy is E n = p n /2m, evolve as free 
particles during the short-time interval between pulses. 
Thus the amplitude to "jump" from the time just after 
the first pulse to just before the second pulse is e~ %Enii '^. 
The final leg of the path shown is another momentum 
jump, caused by the second pulse, from p = p n to p = p m 
and whose amplitude is i m ~ n J m _ n (a). The amplitude to 
proceed from the initial to the final state is the product 
of the amplitudes for the three legs. 



This pathway proceeds from the initial state to the 
final state via the momentum state p n . The system 
can make the transition between these initial and final 
states via any state p n and since these different pathways 
are not detected, the total amplitude for the system to 
jump from the initial to the final state is the coherent 
summation of these individual amplitudes given in Eq. 
(|19l) . This quantum pathways interpretation will enable 
an easy generalization of the two-pulse amplitude to the 
TV-pulse case. 
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D. Comparison with experiment 



i J 



,-i («*)). (34) 



The validity of this theory for two pulses was tested 
in an experiment and reported in Rcf. ([Hj]). In this 
experiment, a BEC consisting of 3 x 10 6 Na atoms con- 
fined in the F = l,mp = —1 ground state by a time- 
averaged orbiting potential (TOP) trap 24[ were released 
and allowed to expand for 1.2 ms as illustrated in the left 
panel of Fig. Next, two 589-nm-wavelength, 100-ns- 
duration, standing-wave, linearly polarized laser pulses 
were applied with a time interval between them which 
varied between 1 and 10 /xs. The intensity of the pulses 
was about 1 W/cm 2 and they were detuned by approxi- 
mately 600 MHz from the 351/2, F = 1 -> 3P 3/2 , F' = 2 
transition. The condensate was then allowed to expand 
for a further 6.2 ms at which time an absorption image 
was taken. This last expansion enabled atoms in non- 
zero momentum states to leave the condensate and the re- 
sulting image is a measurement of the momentum-space 
distribution immediately after the second laser pulse. 

The middle panel in Fig. shows the results of these 
absorption images for varying interval times between the 
pulses. Each row shows a picture of the data for a partic- 
ular interval time while the columns indicate particular 
momentum states. The center column highlighted shows 
atoms in the zero-momentum state. The graph shown on 
the right-hand panel is a comparison of the normalized 
pixel counts of the dots in the zero-momentum states 
(highlighted column in panel (b)) with the theory curve 
|Ao (aoi a o, /3)| 2 plotted as a function of /3 and where ceo 
is the product of the two-photon Rabi frequency and the 
pulse time corresponding to the experimental conditions. 
It is important to note that there are no adjustable pa- 
rameters in this calculation. One can see that there is 
good agreement between theory and experiment which 
in turn lends support for the model presented above. 



IV. TV-PULSE THEORY 

Consider a condensate that is subjected to a sequence 
of TV pulses with arbitrary pulse areas and times. As 
shown in Fig. QJb), we assume that pulse 1 is applied at 
t = t\ and has pulse area a±, pulse 2 at t = i 2 with area 
«2) ■ ■ • , pulse TV at t = t^ with area ajy. If we label 
the momentum state jumped to at the time of pulse k as 
rife, where < k < TV, then a single quantum pathway 
through the entire iV-pulse sequence can be labeled by 
specifying the index of the momentum state the system 
jumps to after each pulse: (rii, 712, . . . , njv-i)- Then it 
is easy to write the amplitude for a particular TV-pulse 
pathway by analogy with the two-pulse case. The result- 
ing amplitude for a single quantum pathway whose start- 
ing momentum state is 2nohk and whose ending state is 
2n^Tik for an TV-pulse sequence labeled in this way is 

Ai N N ] no ( a ,dt,n) = (i^- n °J ni - no {aiY- lE ^ t2 - tl)/n ) 



where we have assumed that the initial momentum state 
is not zero but rather labeled by hq and we have not 
considered any interval following the final pulse at t = t jv. 

The vectors a and St label the and intervals of 

the applied pulse sequence. 



a = (ai,a 2 ,...,a N ) 

St = (Sti = t2 — h, . . . , Stjf-i 



(35) 

tN - t N -i) ■ (36) 



It will be convenient to rescale the intervals in units of 
the Talbot time as was done for two pulses. Thus we 
define fit = S^/Tt for 1 < k < TV — 1. Hence we can 
rewrite Eq. (|3~i)) as 

A^;- n "-iX<x,p) = i n »- n °J ni - no (ai)e- 2vin2lh 

x J„ 2 -„ 1 (a 2 )e- 2 ""^--- 

x e-^-i^J^.^fajv), 

(37) 

To get the full amplitude to jump from the initial mo- 
mentum state, labeled by no, to the final momentum 
state, labeled by njy, we coherently sum over all of the 
single-path amplitudes. This yields the following. 



Ai N \ (a,f3) = Y-Y Ai ni < 

n N ,n a \ ll-J / j / j njv, 



{a, 13) 



m rajv-i 



ni njv-1 



Jn i —no (^1 ) 



■ ■ ■ e" 2 ™ 1 "-^- 1 J„ JV _„ JV _ 1 (a JV ).(38) 

This is the general result for the amplitude to jump from 
an initial momentum state of pi = In^hk to a final mo- 
mentum state of pf = In^Uk due to the application of 
N pulses whose areas are a%, . . . , oln separated by TV — 1 
intervals of durations (expressed in Talbot-time units) 
/?!,..., Pn-x- It holds as long as the Raman-Nath ap- 
proximation is valid for all pulses, i.e., that atoms in 
non-zero-order momentum states do not move an appre- 
ciable distance compared to the condensate size during 
the entire pulse sequence. We can now use this result to 
design sequences of such pulses and intervals to engineer 
a specified momentum-state probability distribution. 



V. ENGINEERING MOMENTUM-STATE 
PROBABILITY DISTRIBUTIONS 

A. Least— squares design 

Designing a pulse sequence to engineer a specified 
probability distribution across the momentum states 
2mhk is straightforward. First, the momentum distribu- 
tion is described by specifying the set of desired probabil- 
ities {p m } for momentum orders 2mhk for all m. We will 







momentum order 



momentum order 



momentum order 



momentum order 



(b) i o. 






momentum order 



momentum order 



momentum order 



momentum order 



FIG. 3: (color online) The figure above shows examples of shaping momentum distributions using only two pulses. The shaping 
parameters are the pulse areas, ct\ and 02, and the interval between them, /3i. All of the plots above show the probabilities 
of atoms being in momentum orders 2mhk where —5 < m < 5. Each plot is annotated with the values of ai, pi, and 0:2 
that produced each momentum distribution, (a) The four plots in the top panel exhibit (from left to right) the distributions 
Dlnag(m) where the two momentum states —2mhk and +2mhk are equally populated and where 1 < m < 4. (b) The four 
plots in the bottom panel depict (again from left to right) distributions Dl 2 J nge (m) where all of the states in the range between 
—2mhk and +2mhk are equally populated and where 1 < m < 4. 



refer to this set of numbers as the momentum probabil- 
ity distribution. This set of numbers must satisfy several 
conditions to be a valid momentum probability distribu- 
tion. Each p m must be a probability and the distribution 
must be normalized so that 



<p m < 1, 



-00 < m < 00, 



and 



]T Pn 

m— — oo 



(39) 



(40) 



Furthermore, since all of the pulses are assumed to be 
standing waves, by symmetry the probability for 2mTik 
must equal that for —2mhk and thus 



Pn 



P- 



—00 < m < 00. 



(41) 



Once the momentum probability distribution is spec- 
ified, the least squares procedure can be carried out by 
defining the TV-pulse, least-squares cost function: 



F 



(N) 



m— — 00 



43 («./*) 



(42) 



and finding the values of the parameters a = ot mm 
and (3 = fj min that produces the global minimum of 

In general, this minimization must be done numeri- 
cally. It is also subject to certain constraints which both 
derive from the Raman-Nath approximation. First, we 



assume that all intervals arc less than or equal to one Tal- 
bot time, or < j3 n < 1 for 1 < n < TV. Second, there 
is a limit on the size of each individual pulse area be- 
cause too large an area will produce population in such 
a high momentum state that there will be appreciable 
motion during the pulse sequence. Thus there is a max- 



imum value, a Tl 



such that < a n < a. 



These 



constraints also limit the total number of pulses that can 
be practically applied while still satisfying all of the con- 
ditions described above. However, as we shall see, many 
important momentum probability distributions can be 
achieved with high fidelity using only two or three pulses. 

We implemented the least-squares minimization pro- 
cedure described above to design two-pulse and three 
pulse sequences to produce momentum distributions in 
two categories. The first category is a distribution where 
all of the population appears in the two momentum states 
where the magnitude is 2nhk (that is, -\-2nhk or —2nhk) 
or equivalently where p n = p—n = 1/2. We will use 

the label, Dmag(n), to stand for the distribution deter- 
mined from the TV-pulse, least-squares procedure when 
the desired distribution is of this type. The other kind of 
distribution we considered has equal populations in the 
range of momentum states between ±2nhk. In this case 
p_ n = ... = p n = l/(2n+ 1). The distribution produced 
by the TV-pulse, least-squares procedure when this type 
of distribution will be labeled by Diange(n). 

These categories of distributions are interesting be- 
cause they act like two- and multiple-beam splitters. 
The ability to transfer condensate population into these 
kinds of momentum distributions might be useful in de- 
signing atom intcrferometry experiments or for quantum 
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FIG. 4: (color online) These plots show examples of momentum distributions shaped by three pulses. The examples shown are 
the same as those shown in Fig. for two pulses. The shaping parameters are the areas of the three pulses, eti, 012, and 03, 
and the intervals between them, /3i and p2- All of the plots above show the probabilities of atoms being in momentum orders 
2mhk where — 10 < m < 10. (a) The four plots in the top panel exhibit (from left to right) the distributions Dmag( m ) where 
the two momentum states —2mhk and +2mhk are equally populated and where 1 < m < 4. (b) The four plots in the bottom 
panel depict (again from left to right) distributions Di^ nge (m) where all of the states in the range between —2mhk and +2mhk 
are equally populated and where 1 < m < 4. 



information processing. In this regard we are particu- 
larly interested in discovering how high a "fidelity" can 
be achieved with just a few pulses. 

Hereafter we will take the measure of the "fidelity," 
that is, how close the actual momentum probability dis- 
tribution is to the specified one, to be the minimum value 
of the least-squares cost function. Thus, 

(f£P) = F[ N S ) (a. m "\ f3 mm ) (43) 

V / min 

will be used to measure the fidelity. 

B. Two— pulse momentum distributions 

For two pulses, the relevant parameters are the dimcn- 
sionlcss pulse areas, a\ and ct2, and the interval between 
the pulses, f3\, measured in units of the Talbot time. 
We obtained least-squares-designed two-pulse sequences 
(which we shall refer to as optimal distributions) for eight 
different specified momentum distributions. These were 
four momentum magnitude distributions whose specified 

(2) (2) 

probabilities arc: D y m ' ag (l) : p_i = pi = 1/2, D\ n ' ag (2) : 

(2) 

P-2 = P2 = 1/2, D y m 'ag(3) '. p- 3 = p 3 = 1/2, and 

(2) 

Dinag(4) P-4 = £>4 = 1/2. As well as four momentum 

(2) 

range distributions: Drange(l) ■ P-i = Pa = Pi = 1/3, 

f2l (2) 
Drange(2) : J5_ 2 = ■ ■ ■ = P2 = 1/5, Drange{l) ■ P-3 = 
(2) 

Pa=P3 = 1/7, DW„,ge(l) ■ P-i = Po = Pi = 1/9. The re- 
sults obtained for all of the optimal distributions in both 
the two-pulse and three-pulse cases are given in Table |T] 
and in Figs. [3] and 0] 



Figure [3^a) shows graphs of two-pulse momentum 
magnitude distributions ordered from left to right along 
the top row. Each graph exhibits the probability for 
atoms to be in each momentum order versus of momen- 
tum orders — 5 < to < 5. The values of the optimal two- 
pulse-sequence parameters and probabilities for each dis- 
tribution arc given in Table |TJ As can be seen from 
the figure, even with two pulses, a remarkably high de- 
gree of fidelity with the desired momentum distribution 
can be obtained with just two pulses. For distribution 
Dmag(l), 43.3% of the population is found in the m = 1 
and m = —1 for a total of 86.6% populating momentum 

(2) 

magnitude 2hk. For distribution Dm' a g(2) we find 42.7% 
of the population in m — 2 and m = —2 states and 35.7% 

(2) 

m m = ±3 (distribution £)^ s (3)). The maximum pop- 
ulation achievable for two pulses degrades to 28.3% for 
to = ±4. We note that this is far better than is possi- 
ble for a single pulse where the probability to populate 
states ±to is J^a). In this case, the maximum values 
arc 33.9% (to = ±1), 23.7% (to = ±2), 18.9% (to = ±3), 
and 16.0% (to = ±4). 

Two pulses do remarkably well in producing equal pop- 
ulations in a range of momentum states. Plots of the op- 
timal distributions determined by least-squares method 
are shown in Fig. EJb) and column two of Table [I] gives 
the values of probabilities achieved for two pulses. One 

(2) 

can see that, for distribution Dra nge (l), the populations 

(2) 

in states to = —1, 0, 1 for distribution D r a nge (\) are the 

approximately the same, p ° ~ p ° pt ^ = 0.3326, to four 
decimal places. This is reasonably close to the specified 

(2) 

value of 0.3333. In distribution _D, Q , lge (2), the popula- 
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TABLE I: The table below presents the optimal values of the parameters and the values of the resulting momentum distribution 
probabilities for two- and three-pulse sequences for the momentum magnitude and momentum range distributions contained 
in Figs. [3] and [4] These values were determined by the least-squares procedure discussed in the text. The leftmost column lists 
the specified momentum distributions in boldface. The next five columns exhibit the parameter values for two-pulse sequences 
while the final seven columns give the three-pulse-sequence results. The two-pulse parameters are the pulse areas, ai and ai, 
and the interval, /3i, expressed in units of the Talbot time. The three-pulse parameters are the areas of the three pulses, ai, 
Qf2, and 03, and the two intervals, fii and pV Columns two and seven, labeled by pi? pt \ give the optimal-distribution values 
of the nonzero specified probabilities. The probability given is listed in the first column in lightface type. Columns six and 
thirteen list the fidelity of the least-squares momentum distributions to the specified distribution by giving the value of the 
least-squares functions F^) and Frg evaluated at the values of the parameters listed in the table. 



Momentum 
Distribution 



Two-pulse results 

~(°P*) ^, rnin omin ~,™i n (J?^)\ 

Tn OL\ El 25 \±LS ) mir} 



Three-pulse results 



P 



(opt) 



_ rain nmin „ rain nmin „ raia / rr>(3) \ 
Ql Pi Q2 P2 Q3 \ F Ls)mi 



P-l =Pl = k 

p (opt) = p (opt) 

P-2 = Pa = | 

(opt) _ (opt) 
P-2 ~ Pi 

P-3 = P3 = | 

(opt) _ (opt) 
P-3 ~ Pz 

P-4 =PA = \ 

P-4, — Pi 
P-l =P =Pl = 

Po 

(opt) (opt) 

P-l =P1 

P-2 = ■ - ■ = Pa = 

Po 

(opt) (opt) 

P-i -Pi 

P-2 ~ P2 
P-3 = ■ ■ ■ = PZ = 

Po 

(opt) (opt) 
P-l =Pl 

P-2 — P2 

(opt) _ (opt) 
P-3 — P3 



-4 = ■ • ■ = P4 = 

Po 

p (opt) = p (opt) 

P-2 — P2 

J^t) — 
P-3 ~ P3 

n (.°P*) — „(<>&) 
P-4 — Pi 



0.4327 

0.4271 

0.3572 

0.2833 

0.3326 
0.3326 

0.1994 
0.1994 
0.1993 

0.1370 
0.1276 
0.1605 
0.1248 

0.0869 
0.1105 
0.1034 
0.1323 
0.0802 



1.715 
2.857 
3.560 
4.230 
1.075 



0.130 
0.320 
0.378 
0.408 
0.153 



0.594 1.5 x 10" 

1.04 1.4 x 10" 

1.429 5.6 x 10" 

1.790 1.2 x 10" 

0.524 3.6 x 10" 



1.861 0.306 0.468 5.6 x 10" 



2.416 0.394 0.836 2.4 x 10" 



0.849 0.218 3.895 4.8 x 10" 



0.4766 

0.4877 

0.4179 

0.3984 

0.3328 
0.3328 

0.1998 
0.1998 
0.1998 

0.1402 
0.1410 
0.1423 
0.1405 

0.1087 
0.1108 
0.1087 
0.1120 
0.1073 



3.737 0.481 3.402 0.453 1.782 1.9 x lO"- 3 

2.733 0.342 0.585 0.337 0.882 3.9 x 10~ 4 

1.297 0.586 7.347 0.404 2.310 1.7 x 10~ 2 

3.158 0.376 2.936 0.287 1.047 2.5 x 10~ 2 

0.551 0.565 1.676 0.856 0.547 1.4 x 10~ 6 

1.730 0.058 0.206 0.252 0.491 5.0 x 10~ 7 

2.205 0.451 0.739 0.260 2.112 6.7 x 10~ 5 



2.436 0.267 1.627 0.402 0.949 1.3 x 10" 



tions in states m = — 2, — 1, 0, 1, 2 the distribution across 
these states vary between 0.1993 and 0.1994 which com- 
pares well with the specified value of 0.2. For distribu- 

(2) 

tion, Drange{3), while the states m = —3, —2, —1, 0, 1, 2, 3 
contain more than 96% of the total population, the prob- 
abilities vary between 0.1248 and 0.1605 producing a 
relatively large variance around the specified value of 

1/7=0.1429. The degradation of the fidelity can espe- 

(2) 

cially be seen in distribution Drange (4) both in the value 
of (F^g) m i n and in the deviation of the achievable prob- 
abilities from the specified probability as shown in Table 
HI For this distribution, the probabilites vary between 
0.0802 and 0.1323 and deviate significantly from the spec- 
ified value of 1/9=0.1111. 



Upon closer inspection, we found that the optimal two- 
pulse sequences had the common characteristic that there 
was a single dominant pathway to each of the prescribed 
final momentum states. Thus, for a fixed final momen- 
tum state, this dominant pathway consisted of a direct 
jump from the zero-momentum original condensate to 
the final momentum at the first pulse followed by evolu- 
tion between the pulses and no jump in momentum at the 
second pulse. In this single-dominant-pathway picture, 
the total amplitude is somewhat insensitive to the time 
between pulses since the coherent sum can be (roughly) 
approximated with a single term. 
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C. Three— pulse momentum distributions 

For three pulses, the relevant parameters are the pulse 
areas, a±, and 03, and the intervals between the 
pulses, /?i and (3 2 - We obtained optimal three-pulse se- 
quences for the same set of eight specified momentum 
distributions as for two pulses. As can be seen in Fig. [4l 
the fidelities achievable with three pulses is better than 
with two pulses and is far superior to the single pulse 
case. 

Figure EJa) (top row) shows three-pulse momentum 
magnitude distributions D y m ' ag (l), D y m ' ag (2), D y m ' ag (3), 

and -Dma g (4) respectively from left to right. While the 
bottom row shows three-pulse momentum range distri- 
butions, Dr 3 ange(l), -D™L,e(2), £>ra„ ff e(3), and D^ange^) 

again from left to right. These plots show that the three- 
pulse optimal distributions faithfully reproduce the spec- 
ified distributions better than the two-pulse versions. 
Comparisons of the fidelities for the two-pulse case versus 
the three-pulse given in Table Q] shows that three pulse 
does a better job at reproducing the specified distribution 
in every case. 

For some of the distributions there is little difference 
between two pulses and three pulses because of the high 

(2) 

fidelity of the two-pulse case. One example is Dra nge (l) 
versus Dr a ' nge (l). However, three pulses is clearly better 
for the -Dran 9 e(3) and Drange (4) distributions. 

VI. CONCLUSION 

In this paper we have shown that it is possible 
to create clouds of coherent atoms with momentum- 
space distributions important for applications in atom- 
interferometry with a sequence of only two or three 
standing-wave laser pulses applied to a BEC. We de- 
rived the momentum distribution for atoms in a BEC 
after TV short-time, standing-wave laser pulses were ap- 
plied where the pulse areas and time intervals between 
the pulses were variable. This distribution was a gener- 
alization of the expression for two pulses whose valid- 
ity has been verified experimentally [ll|. We further 
described a method for designing pulse sequences that 
produce a specified momentum distribution of the con- 
densate atoms. We found that two kinds of distribu- 
tions that have important applications as beam splitters 
can be produced with high fidelity with two or three 
pulses. We also found that the optimal two-pulse se- 
quences obtained could be understood in terms of a 
single-dominant-pathway picture. 

The ability to produce coherent atom with engineered 
momentum-space distributions can now become a new 
tool for the design of new atom interferometer schemes. 
Methods for initial momentum-state selection or for pro- 
ducing multiple-beam splitters can now be designed. For 
example, one could imagine a Bose-Einstein condensate 
created and confined on an atom chip to which could be 



applied a sequence of pulses such as Dr an ge(2) so that the 
cloud is split into five equal parts. If such a condensate 
were confined by a harmonic trap potential, these parts 
would eventually all come back together at once where 
they could be split again. This would produce multiple 
interference patterns that reflect the different phase evo- 
lutions along the different pathways. Such multi-particle 
interferometers could, in principle, implement quantum 
computations or be used for precision navigation appli- 
cations, gradiometry, or fundamental studies. 

Appendix A 

This appendix presents the derivation of Eq. ([9]) which 
expresses the condensate wave function just after the ap- 
plication of a short-time, standing-wave laser pulse in 
terms of the wave function just before the pulse. 

We begin by transforming away the internal energies 
of the atom: 

^k(r„,i) = e-' lEkt / h 4> k {v a ,t). k = g,e (Al) 
Under this transformation, Eqs. ([7]) and © become 

+ hn cos (k x • r„) / (t) e.- 4A Ve (r a , t) 

+ 5 A(j0 3 (r a ,i)| 2 + |0 e (r a ,t)| 2 ) (f> g (r a ,t), 

(A2) 

and 

+ MoCOs(k L -r a )/(t)e iAt s (ra,O 

+ , 9 A^(|0 g (r a ,t)| 2 + |0 e (r a ,t)| 2 ) (f> B (r a ,t), 

(A3) 

where we have made the rotating-wave approximation 
(RWA). The RWA consists of neglecting the exponen- 
tials that oscillate as exp [±i (ujq + wl) t] relative to 
exp(±iAi) where Hujq = E e — E g . This approximation 
holds because wi is typically six orders of magnitude 
larger than A for optical frequencies when the light is 
close to resonance. 

Since we assume fio <C A, there is never very much 
population in the upper state. Thus we can neglect the 
nonlinear term and kinetic plus trap potential energy 
terms in Eq. (|A3[) and can neglect <j> e in the nonlinear 
term of Eq. (|A"2j) . This gives 

ih^f- « ^ 9 V<,(ra,t)+gA|0 g (r Q ,t)| 2 g (r a ,t) 
+ hQ cos (k L • r a ) / (t) e~ iAt 4> e (r„, t) 

(A4) 
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ih-g « ftfi cos(k L T a )/(i)e iA V ff (r„,t) (A5) 

Finally, note that the first two terms on the right-hand- 
side of Eq. (|A4[) approximately satisfy the static GP 
equation and thus together equal fJt(f)g. Replacing those 
two terms gives: 



+ hn o cos(k L -r a )f(t)e~ lAt e (r a ,t) 

(A6) 

This term can be transformed away but it is unecessary 
since we have assumed that fj,St/h <C 1. Neglecting this 
term finally gives 

ih^f « ftfi cos(k L T a )/(t)e- <A Ve(ra,t)(A7) 



ift-^ « M} cos(k L -r a )/(t) e * A V 9 (r a ,t) (A8) 

Under the "sudden approximation," these equations can 
be easily solved. That is, we assume that the turn-on 
of the pulse is fast enough that the atom remains in its 
initial state until the pulse is fully on. 

In order to express the solution in terms of the area of 
a single pulse, we transform the solutions as 

(d> g (T a ,t)\ = (e-**/* \($ g (v a ,t)\ 
\4>e{v a ,t)) I e' A */ 2 ) Ue(r 0! t) ) 



This yields the following equations for 



(A9) 



l% ^tr = ^ A ^(ro,*) + V(r a ,t)^(r OI t)GA.10) 



where 



Defining 



^(r ,t) = ftfio/(t)co8(k £ .r ). (All) 



and 



U{r a ,t) = 



I <Pe (r a ,t) 



■±hA V(r a ,t) 



V(r a ,t) iftA J 
we may formally express Eqs. (|A10|) as 
<9<l 



(A12) 



(A13) 



(A14) 



Since H (r a ,i) is constant during the pulse we obtain a 
relationship between $ before and after the pulse 

$ (r a , to + St/2) = e - iM5t $ (r„ t - St/2) , (A15) 

where 



M _ | 2 " ^ cos(k L • r a ) 

I fi cos(k L • r a ) ±A 



(A16) 



This matrix can be easily exponentiated using its eigen- 
values, 



A 



± 



(iA) 2 + ft 2 cos 2 ( kL . ra ) 



1/2 



= ±A, 



(A17) 



and eigenvectors 
l+A> = 



sin (0/2) 
cos (6»/2) 



A) 



cos (6>/2) 
-sin (0/2) 



where is defined by 

O cos(k L ■ r ) 



sin (0) 



A 



cos (0) 



A ' 



The diagonalization matrix, U* , where 



Mn = UMU^ = 



A 

-A 



is given by 



[/t = / sin (9/2) cos {9/2) 
\ cos (6»/2) -sin (9/2) 

The exponentiated matrix is, therefore 



(A18) 
(A19) 

(A20) 
(A21) 



-iMSt 



g 2 g— i\St 



c 2 e iXSt 



iXSt 



iXSt 



-sc (e iXSt — e - ^ 5 *) 



-sc (e 

— iXSt _|_ ^2giA(5# 



(A22) 

where s = sin (9/ '2) and c = cos (9/2). 

Thus we can write a relationship between </> g before 
and after the pulse. From Eqs. (|A15[) and (|A22|) we have 



(r„,*+)= (s 2 e- iX5t + " 2 - iXSt 



sc (e iX5t -e~ iX5t 



) 4>_ g (r .*-) 
) e (r Q ,t_) 



(A23) 



where t± = to ± St/2. 

We can simplify the above expression by invoking the 
approximation that the detuning from resonance is much 
larger that the single-photon Rabi frequency, !! « A. 
In this case, we have 

A « ^A + cos 2 (k L • r a ) 

= iA + i^](l + cos(2k L .r Q )) (A24) 
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Furthermore, 



to set s~0 and c ~ 1. Hence 



cos(0) = V 



iA) 2 + fi2 C os2 (k L -r a )) 1/2 
l-^(T|) 2 cos 2 (k L .r Q ). (A25) 



The factors s 2 and c 2 are, to second order in CIq/A, ap- 
proximately 



s 2 = sin 2 (<9/2) 
1 



-(i-c Os( 0)) = i(i-i£ 



2 _ „ 2 



(A26) 



c 2 = cos 2 (6»/2) 



2 (l + c 0sW ) = |ll + ^ 



A 



' "° * cos 2 (k L • r„) . (A27) 



c 1 a 1 



Given the above it should be an excellent approximation 



(r„,t+) ra e i(iA+n 2 )«+m 3 « 003(2^.^)^ ( fojt _) ) 

(A28) 



where 



f2 2 
2A 



(A29) 



is the two-photon Rabi frequency. We can express the 
above relationship in terms of the original condensate 
wave function as 



4> g (r ai t + ) w e 

x e m2St cos{2kL "») fl (r„ , *_) . (A30) 

This is the effect of a single standing-wave pulse on the 
condensate wave function. 
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